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ABSTRACT 

The continuing difficulty of achieving a reliable explosion in simulations of core-collapse supernovae, 
especially for more massive stars, has led to speculation concerning the observable transients that 
might be produced if such a supernova fails. Even if a prompt outgoing shock fails to form in a 
collapsing presupernova star, one must still consider the hydrodynamic response of the star to the 
abrupt lo s s of m ass via neutrinos as the core forms a protoneutron star. Following a suggestion by 
Nadezhin (1980), we calculate the hydrodynamical responses of typical supernova progenitor stars 
to the rapid loss of approximately 0.2 to 0.5 Mq of gravitational mass from their centers. In a red 
supergiant star, a very weak supernova with total kinetic energy ~ 10 47 erg results. The binding 
energy of a large fraction of the hydrogen envelope before the explosion is of the same order and, 
depending upon assumptions regarding the neutrino loss rates, most of it is ejected. Ejection speeds 
are ~ 100 km s _1 and luminosities ~ 10 39 erg s _1 are maintained for about a year. A significant part 
of the energy comes from the recombination of hydrogen. The color of the explosion is extremely red 
and the events bear some similarity to "luminous red novae," but have much lower speeds. 

Subject headings: black holes - supernovae: general, stars:neutron, neutrinos, shock waves 



1. INTRODUCTION 

Black holes are expected to fo rm in a significant frac- 
tion of massive star deaths (e.g.. lO'Connor fe Qttll20lTI : 
lUgliano et al.1120121) and it seems possible that in at least 
some of these cases the black hole will form without gen- 
erating an outgoing shock. What observable signal would 
accompany such an eve nt? Would they all be "un- novae" 
(jKochanek et al.l 120081 ) . stars that just disappear from 
view, or will some sort of transient display announce the 
birth of any black hole? 

If the progenitor star was rotating, then a variety of en- 
ergetic transients might be possible depending on the dis- 
tribution of angular momentum, ranging from common 
gamma-ray bursts by the collapsar mechanism (jWooslevi 



the total energy emitted should be close to the binding 
energy of a neutron star with its maximum stable gravi- 
tational mass, which we allow to vary in the range 2.0 - 
2.5 M Q in our model. 

As frequently noted, this huge energy loss also results 
in a decrease in the gravitational mass of the compact ob - 
ject (jLattimer fc Yahillll989HLattimer fc Prakashll2001l) . 
This mass loss is of order 
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(jQuataert fe ShiodeH2012t IWooslev fe Hegerll2011l) tran- 
sients. But what if the star were rotating very slowly or 
not at all? Pulsation s (IWooslev fc Hegerf 2007) or acous- 
tic energy transport (jQuataert fc Shiodell2012j ) may lead 
to envelope ejection, but these either occur in stars that 
are uncommonly massive or depend on details of gravity 
wave propagation that are still to be worked out . 

Here we follow on a suggestion by INadezhinl (|l980) 
which provides a simple mechanism for ejecting at least 
some mass in supernova progenitors with very weakly 
bound envelopes, i.e., red supergiants. In all but ex- 
tremely rare super-massive stars , iron core collapse lead s 
to protoneutron star formation ([O'Connor fc Ottl feOll). 
The protoneutron star must then emit a considerable 
fraction of its mass-energy as neutrinos before either set- 
tling down as a cold degenerate object or contracting 
inside its event horizon and becoming a black hole. This 
is true even if the protoneutron star continues to accrete 
substantial mass while contracting. In fact, in that case, 
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where Mq is the gravitational mass of the cold neutron 
star. This implies a mass loss of approximately 0.2 - 
0.5 Mq over a period of several seconds, far shorter than 
the sound crossing time for the helium and heavy element 
core. This "mass" streams out through the presupernova 
star at the speed of light without interacting appreciably. 
To the remaining star it simply appears that the gravi- 
tational potential of the core has abruptly decreased. In 
response, the star begins to expand. In this paper we 
simulate this mass loss and the resulting expansion to 
determine whether it is capable of creating a shock with 
sufficient energy to reach the outer layers of the star and, 
if it does, to eject mass and form a visible transient. It is 
important to note that this expansion precedes the loss 
in support pressure caused by the collapse of the core. 
That rarefaction passes out through the star on a slower, 
free-fall time scale. 

2. COMPUTATIONAL PROCEDURE 

We began with two presupernova models calculated 
using the KEPLER stellar evolution code. KEPLER 
is a Lagrangian ID implicit hydrodynamics code with 
the appropriate nuclear physics, mass loss, opacities, 
and eq uations of state for studying massive st ar evo- 
lution d^^reFn] EOSS Imosi^SrS] iOOit). It is 
not Courant-limited, allowing it to follow a star from 
its formation to the collapse of its iron core. How- 
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ever, because KEPLER uses a Lagrangian grid, it can- 
not easily implement an absorbing inner boundary con- 
dition. Consequently, it cannot accurately follow the 
infall of the star once a compact object has formed at 
its center. Once the presupernova models calculated 
by KEPLER reached supersonic collapse speeds in their 
cores, they were thus mapped into CASTRC0 for fur- 
ther study. CASTRO is a multidimensional Eulerian 
radiation-hydrodynamics code with adaptive mesh re- 
finement, stellar equati ons of state, nu clear reaction net- 
works and self-gravity (|Almgrenll2010t ). The calculations 
were done in one dimension with reactions turned off, and 
the central object modeled as a point of variable mass. 
An Eulerian grid with constant mesh and spherical ge- 
ometry was employed. After the collapse had proceeded 
long enough in CASTRO to allow shock propagation to 
the base of the hydrogen envelope, the results were then 
mapped back into KEPLER again and the final hydro- 
dynamics and a light curve were calculated. 

2.1. Central Object Modeling 

Hydrodynamics and radiation transport in the region 
immediately surrounding a protoneutron star or black 
hole is complex and difficult to model. A realistic sim- 
ulation of the neutrino transport alone is well beyond 
the scope of this paper and the outcome would depend 
upon the chosen neutron star equation of state and the 
dimensionality of the calculation. Since we are only in- 
terested in the temporal evolution of the central gravita- 
tional potential, however, reasonably good results can be 
obtained from parametrized calculations with a generic, 
qualitative description of the neutrino energy loss. The 
inner boundary of our simulation is placed, in CASTRO, 
at the outer edge of the pre-collapse iron core. Matter 
interior to this boundary is treated as a point having 
variable gravitational mass. Initially this mass is equal 
to the baryonic mass of the interior matter, i.e., the iron 
core mass. As time passes that gravitational mass de- 
creases due to neutrino emission and, if no mass were 
added, it would eventually decrease to the known cold 
neutron star value for a given ba ryonic mass and equa- 
tion of state (|Prakash et al.l [l997). However, appreciable 
mass does accrete, leading in the more interesting cases 
to the formation of a black hole. During the accretion 
phase neutrinos continue to carry away effective mass. 
Since we neglect rotation, we assume that, after the black 
hole forms, matter and energy flow into the black hole 
without any more emission. 

An upper limit to the mass lost is given by the gravi- 
tational binding energy of the maximum stable neutron 
star mass, the "Tolman-Oppenheim er-Volkoff (TOV) 
limit" (|Oppenheimer fc Volkofl 1939). Current best lim- 
its place this param eter in the range 2.0 to 2.5 M^of 
gravi tational mass (|Akmal et al.l 119981 : iDemorest et al.l 
2010). The characteristic time scale for this loss, which 
is assumed to occur exponentially, is given by the cool- 
ing t imescale parameter 77.., tak en here to be 3 seconds 
(e.g. iBurrows fc Lattimerlfl9 87). For this simple case, 
the mass loss rate is 



This equation gives an upper bound on the mass loss and 
therefore on the strength of the transient produced. This 
model is referred to as the "maximum mass loss" case. 

Somewhat more realistically, we can take into account 
the time for the critical mass to accrete and follow the 
neutrino energy losses that occur during that time. For 
each bit of mass that accretes, some fraction of its grav- 
itational binding energy is radiated away promptly, but 
more is emitted with a delay since the neutrinos must dif- 
fuse out. Once the black hole has formed, the neutrino- 
emitting material falls within the event horizon much 
faster than the neutrino diffusion timescale, and any en- 
ergy that has not been rel eased by the time of collapse 
ends up in the black hole ()Prakash et al.lH997l : iBurrowsl 
1988). The protoneutron star is also extremely hot and 
thus behaves in a slightly different manner from cold neu- 
tron stars. We represent the protoneutron star as a low- 
entropy core surrounded by a hot envelope and track 
three masses: Men, Mb, and Mth- 

Mb is the baryonic mass, i.e. all mass that either was 
inside the inner boundary at the time of collapse or sub- 
sequently accretes through the inner boundary. Mch is 
the effective gravitational mass of the hot neutron star, 
which is not, in general, equivalent to M a , the gravita- 
tional mass of a cold neutron star with the same baryon 
content. This is the mass that defines the gravitational 
potential in the outer part of the star. M t h accounts for 
the difference in mass-energy stored in this hot neutron 
star as opposed to a cold one, due to both its higher 
internal energy and inflated radius. This extra internal 
energy diffuses away on the cooling timescale r c . The 
gravitational mass Mch of the hot neutron star is then: 
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M Gh = M + M th (3) 
= M B - BE C + M th (4) 

where BE C is the cold neutron star binding energy. As 
matter accretes a fraction e of the subsequent change 
in binding energy is assumed to be trapped and added 
to M t h while the remaining fraction 1 - e is radiated 
promptly as neutrinos. The trapped energy diffuses out 
on a cooling timescale r c . The evolution of M t h is thus 
given by 

sV Mth 1 dBE c ■ 

Mth = + e -jr-j- M B 5 

t c dM B 

and the initial condition M t h = BE C . The derivative 
d BE c I dMB is evaluated from the binding energy equa- 
tion, Eq. [TJ Once all Mth has diffused away, i.e. the 
neutron star has cooled, Eq. 0] reduces to the standard 
cold neutron star equation Mq = Mb — BE C . Assuming 
rapid virialization of the accreted material implies that 
the parameter e should not exceed ~^0.5. As the cooling 
timescale at the surface is likely less than the timescale 
for the full PNS, e is likely to be much lower. In the 
limit that the cooling timescale were much shorter than 
the accretion timescale, e would go effectiv ely t o zero. 
For our models we adopt e = 0.1. In Section [3~3l we test 
the effects of varying this parameter and find them to be 
small. 

The emission halts when the neutron star has accreted 
past the TOV maximum mass limit. However this limit- 
ing mass must be adjusted for the neutron star's evolv- 
ing heat content. Therefore we wait until the cold core 
of the neutron star (Mch — Mth) nas accreted past the 
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Fig. 1. — Core mass growth in the 15 Mq model RSG15 with an 
assumed TOV mass limit of 2.5 Mq, demonstrating the "full loss" 
model. The black curve (solid) shows Mq^, the gravitational mass. 
Blue (dash-triple dot) is Mb, the baryonic mass. Red (dash-dot) 
is the thermal mass M t ^. Green (dashed) shows the cumulative 
lost mass from the simulation. 

TOV limit before we assume the core has become a black 
hole and shut off the neutrino emission. From this point 
on the central object behaves as a purely gravitational 
point mass that absorbs all of Mb ■ This is referred to as 
the "full loss" model. 

Figure[T]shows an example of the full model for RSG15, 
TOV = 2.5 M Q . The blue curve shows M B , the total 
baryonic mass that has entered the core region (prec- 
ollapse core plus accretion); it is higher than the black 
curve showing Mchi the effective gravitational mass of 
the core. The red curve shows M t h, which decays ex- 
ponentially as the protoneutron star cools. The green 
curve shows the cumulative mass loss in the simulation. 
Even though the core in this case lives for almost 24 s as 
a neutron star, most of the mass loss occurs during the 
first 5 seconds. 

2.2. Choice and Structure of Stellar Models 

The two presupernova models used for this paper were 
red supergiants with a ZAMS mass of 15 M Q (RSG15), 
shown in Fig. 03 and 25 M Q (RSG25), shown in Fig. 
[31 both of solar metallicity . The y are taken from the 
survey of iWooslev fc Hegen (|2007l) . Both stars shed sub- 
stantial amounts of mass before reaching the end of their 
lives. The 15 Mq star represents a more common su- 
pernova progenitor, while the 25 Mq model has a more 
compact structure and may be more likely to make a 
black hole promptly. The observed upper limit on the 
progenitor masses of standar d core-collapse supernovae 
found to date is 18 M (e.g. iSmarttl 120091 ). suggesting 
that higher mass progenitors are more likely to collapse 
in a non-traditional way, or that these progenitors are 
preferentially obscured. Model RSG15 has a helium core 
of 4.27 M that extends to 3.568 x 10 10 cm. RSG25 has 
a helium core of 8.20 M Q that extends to 1.807 x 10 10 
cm. Both red supergiants have an extended, tenuously- 
bound hydrogen envelope extending from the helium core 
out to ~ 5 x 10 13 cm that is easily ejected. The net bind- 
ing energy of the envelope external to what is nominally 
the helium core (where the hydrogen mass fraction de- 
clines to 1%) is 1.1 x 10 48 erg, but a short distance out at 
4.48 M this declines to 10 47 erg, a value that character- 



1.000 



0.100 



0.010 



0.001 




H ~ 

He4 

C12 - 

N14 

016 _ 

Ne20 

Mg24 - 



Ca40 — 



10° 



10' 1 



10" 



r (cm) 



Fig. 2. — Composition of the RSG15 supernova progenitor star. 
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Fig. 3. — Composition of the RSG25 supernova progenitor star. 

izes most of the hydrogen envelope. These values include 
internal energy but not the energy potentially available 
from recombination. For the 25 Mq model, the net bind- 
ing energy external to the helium core is 6.4 x 10 erg, 
but this declines to 10 47 erg when the interior mass in- 
creases from 8.20 Mq to 8.36 M Q and a still smaller value 
characterizes most of the envelope. We note that the re- 
combination of 10 Mq of hydrogen (13.6 eV per atom) 
would release 2.6 x 10 47 erg. 

Some additional parameters of both models can be 
found in Table 1. Also g iven is the compactnes s pa- 
rameter, £ 2 .5, as defined bv lO'Connor fc OttJ (|2011l ). but 
computed at the time when the star is moved to CAS- 
TRO rather than at the time of core bounce. A star with 
a higher compactness parameter has a denser region sur- 
rounding its core; for our purposes, this means it will 
accrete to the TOV limit faster and therefore potentially 
lose less mass as neutrinos. 

Because of the extended size of the red supergiant pro- 
genitors, we do not carry the entire star on our simulation 
grid in CASTRO. Instead we model only the helium core 
and base of the hydrogen envelope. 

3. HYDRODYNAMIC RESPONSE 

3.1. Models Without Mass Loss 

Without an outgoing shock or mass loss from neutri- 
nos, the inner layers of the star should collapse directly 
into the black hole. This scenario provides an excellent 
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TABLE 1 
Stellar Model Parameters 



Model Name 


Final Mass (M ) 


He Core Mass (M ) 


He Core Radius (cm) 


Compactness §2.5 


RSG15 


12.79 


4.27 


3.568X10 10 


0.18 


RSG25 


15.84 


8.20 


1.807X10 10 


0.33 



400 



200 - 




Fig. 4. — Velocity curves showing collapse of the RSG15 stellar 
model in the prompt black hole formation case, i.e. without any 
mass loss. Positive velocities are outwards, negative are inwards. 
Curves are purple for early times, shading to red for late, approxi- 
mately 15 s spacing. With no mass decrement and no core bounce 
shock to provide outward velocity, the star simply falls into the 
black hole. Total time shown: 709 s. 

check of the fidelity of our simulation. As Figure [4] shows, 
our model accurately reproduces this behavior. Dark 
(purple) colors on the plot indicate early times, shading 
to light (red) colors at late times. Initially the bulk of 
the star is in hydrostatic equilibrium (zero velocity) with 
a small portion near the core showing high infall veloc- 
ities. As the collapse continues, more and more of the 
stellar material acquires negative velocities. Without an 
outgoing shock to counter this motion, it eventually falls 
into the core. If the simulation is run for long enough, all 
mass will disappear from the grid. Physically, this sce- 
nario corresponds to prompt black hole formation, when 
the core collapses without passing through an interme- 
diate protoneutron star stage. If the outer layers of the 
star do not have sufficient angular momentum to form a 
disk when they reach the black hole, the entire star will 
disappear without producing a transien t - the unnova 
case as defined by iKochanek et al.l (|2008|) . 

Figure [5] shows the growth of the simulated point mass 
in RSG15. The black curves show the model with no neu- 
trino mass loss; in this case the point mass corresponds 
to Mb- The dashed line represents the growth of the 
point mass while the dot-dashed line shows the mass on 
the simulated grid. Over time this point mass grows as 
the grid mass declines, with most of the change occurring 
in the first 15 seconds. The solid line shows the sum of 
the point mass and grid mass, i.e. the total mass rep- 
resented in the simulation; in the no-mass-loss case it is 
constant throughout. This confirms that our simulation 
is reproducing the collapse accurately, without spurious 
shocks or unphysical mass loss. 

3.2. Models With Maximum Mass Loss 
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Fig. 5. — Various core masses as a function of time for Model 
RSG15 and a TOV limit of 2.5 M . The central mass M Gh 
(dashed) grows while the mass on the simulation grid (dot-dashed) 
drops. The solid line represents the total mass in the center and on 
the grid combined. Black curves show the no neutrino loss (con- 
stant mass) case. Blue curves show the maximum-loss model, while 
the red curves show the full neutrino loss model. 

The set of blue curves in Figure [5] shows the evolu- 
tion using the simplified maximum mass loss model. In 
this case, the point mass corresponds to Mq with losses 
defined by Eq. [2 in which the core loses the total bind- 
ing energy appropriate to a TOV-limit neutron star on 
a time scale t c , regardless of the amount of mass flow- 
ing in from the collapsing star. This model therefore 
provides an upper bound on possible transients, as the 
core cannot lose more mass than the binding energy of 
the largest possible neutron star. The amount of mass 
lost in each case is listed in Table 2. Following the blue 
curves in Figure O as the collapse begins the point mass 
growth (dashed) is noticeably suppressed for the first 5 
seconds as neutrino losses balance accreted mass. The 
overall mass in the simulation (solid) drops accordingly, 
then becomes constant as mass loss ceases. 

This mass decrement was sufficient to produce an out- 
going shock in the inner layers of the 15 M Q presupernova 
star. The shock's evolution is seen in Figure |6] (purple 
at early times, shading to red at late times). The shock 
grows in speed as it leaves the helium core, and succeeds 
in reaching the base of the hydrogen envelope. Of inter- 
est is the fact that the shock strength varied noticeably 
with the choice of neutron star upper mass limit. The 
approximate shock strengths at 10 11 cm for our six dif- 
ferent choices of TOV limit are listed in Table 2. 

3.3. Models With Realistic Neutrino Mass Loss 

The red curves in Figure [5] show the mass evolution 
for the full model for neutrino losses described in Section 
2.1. In this case the dashed line showing the point mass 
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Fig. 6. — Collapse of the RSG15 model in the maximum neutrino 
mass loss case, showing a shock forming due to the effective core 
mass decrement. Positive velocities are outwards, negative are in- 
wards. Curves are purple for early times, shading to red for late. 
This model has the TOV limit set to 2.5, resulting in a mass decre- 
ment of 0.525 Mq. The shock propagates out of the helium core 
(r = 3.568xl0 10 cm). Time shock reaches 10 10 cm: 38 s. Time to 
end of helium core: 196 s. Time to 10 11 cm: 577 s. 

corresponds to Man as given by Eq. |U This model takes 
into account the thermal mass loss and ties the cessation 
of neutrino losses to the amount of material that has been 
accreted by the core, rather than switching it off after a 
predetermined timescale as in the upper bound model. 
As can be seen in Figure [SJ this model (red) loses mass 
over a longer timescale than the maximum loss model 
(blue) , continuing until the point mass reaches the TOV 
limit, in this case 2.5 Mq, after which the total mass 
becomes constant. We expect equal or less mass loss in 
this case as compared to the maximum loss model. In 
the case where the TOV limit is high enough that the 
neutron star lives for longer than the cooling timescale, 
the core loses close to the maximum possible amount of 
mass; in the case where it does not, however, mass loss 
is suppressed as neutrinos that would have been emitted 
instead end up inside the black hole. The amount of mass 
lost in each case is listed in Table 3. 

Though the overall mass decrement in the full model 
cases is lower than in the maximum loss case, it is still 
sufficient to produce an outgoing shock. Figure [7] shows 
the shock evolution for RSG15, TOV = 2.5 M©. The 
approximate shock strengths at 10 11 cm for our six dif- 
ferent choices of TOV limit are listed in Table 3. The six 
shocks created in RSG15 are shown in Figure [5J 

We also tested variations in the parameter e, which 
controls the fraction of binding energy trapped as ther- 
mal mass. Changes in e have a small but real effect on 
the total mass loss, depending on the amount of accreted 
mass. The more mass accreted, the more important e will 
be. As higher TOV limit models tend to accrete longer, 
e has a higher impact here. A lower epsilon leads to a 
higher mass loss as less of the binding energy is tem- 
porarily trapped as thermal mass. We tested the range 
e < 0.5, identified as the physically reasonable range of 
this parameter. For the case where the TOV limit is 2.5 
M©, the most sensitive, a change of 0.05 in e in RSG25 
resulted in approximately a 0.011 M change in the over- 
all mass loss. In the extreme case e = 0.5, this will make 



Fig. 7. — Collapse of the RSG15 model in the fully-modeled neu- 
trino mass loss case, showing a shock forming due to the effective 
core mass decrement. Positive velocities are outwards, negative 
are inwards. Curves are purple for early times, shading to red for 
late. This model has the TOV limit set to 2.5, resulting in a mass 
decrement of 0.523 Mq. The shock is smaller in strength than the 
maximum-loss case and reaches the edge of the simulation with a 
lower velocity. Time shock reaches 10 cm: 40 s. Time to end of 
helium core: 207 s. Time to 10 11 cm: 620 s. 
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Fig. 8.— RSG15 shocks at the limit of the CASTRO simulated 
domain for six different choices of TOV limit, full neutrino loss 
model. 

a TOV = 2.5 M model look like a TOV - 2.35 M 
model. For lower TOV limits, the effect of varying e 
was smaller. Thus, although a change in e will affect the 
mass decrement and by extension the shock strength, the 
results are robust. 

3.4. Effects of Stellar Structure 

The kinetic energy of the shock that reaches the hy- 
drogen envelope is strongly influenced by the size of 
the presupernova core and the overlying stellar struc- 
ture through which it must travel, as can be seen by the 
differences in kinetic energy between RSG25 and RSG15 
in the maximum mass loss model. This model loses the 
same amount of mass in the same time regardless of stel- 
lar structure, but the final shock in RSG25 models is 
significantly weaker than in RSG15, having propagated 
through a much heavier carbon-oxygen and helium core 
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TABLE 2 
Maximum Mass Loss Models 



Stellar Model TOV (M ) Mass Lost (M ) KE a (ergs) Shock Strength b (km/s) 



RSG15 2.0 0.336 1.875xl0 47 902 

2.1 0.370 2.599X10 47 985 

2.2 0.407 3.572X10 47 1070 

2.3 0.444 4.855X10 47 1158 

2.4 0.484 6.554X10 47 1249 

2.5 0.525 8.719X10 47 1341 



RSG25 2.0 
2.1 
2.2 
2.3 
2.4 
2.5 

a At base of hydrogen envelope 
b At r = 10 11 cm 



TABLE 3 
Full Mass Loss Models 



Stellar Model TOV (M ) Mass Lost (M ) KE a (ergs) Shock Strength b (km/s) 



RSG15 2.0 0.277 1.287xl0 47 814 

2.1 0.331 2.059X10 47 926 

2.2 0.382 2.953X10 47 1019 

2.3 0.430 3.911X10 47 1094 

2.4 0.477 4.896X10 47 1157 

2.5 0.523 5.779X10 47 1204 



RSG25 2.0 
2.1 
2.2 
2.3 
2.4 
2.5 

a At base of hydrogen envelope 
b At r = 10 11 cm 

compactness) and a more massive iron core of 1.83 M , 
compared to 1.63 M . As it therefore accretes faster 
than RSG15 and is already closer to the TOV limit, the 
core spends significantly less time in the protoneutron 
star state; consequently the same choice of parameters in 
the larger star causes less mass loss. The full mass loss 
model in RSG25 also shows systematically lower mass 
decrements than the maximum mass loss model in all 
cases, indicating that even with a high TOV limit and a 
relatively long-lived protoneutron star, not all the bind- 
ing energy is being emitted before collapse to a black 
hole. The six shocks produced in RSG25 are shown in 
Figure El 

4. TRANSIENTS PRODUCED 

In the case of prompt black hole formation in a star 
without a high-J outer layer, there is no shock and no 
visible transient - the star simply disappears. Here we 
consider further the case of protoneutron star formation 
and delayed collapse, in which an outgoing shock has 
formed. Will a detectable transient be produced, inter- 
mediate between a complete d isappearance and a tra - 
ditional explosion, as raised by iKochanek et ail ((2008)? 
To evaluate this question we took models from CASTRO 
where the shock had reached 1 x 10 11 cm (the limit of 
the simulation) and mapped them back into KEPLER, 



0.336 6.537X10 46 723 

0.370 1.002X10 47 820 

0.407 1.483X10 47 919 

0.444 2.139X10 47 1025 

0.484 3.011 xlO 47 1134 

0.525 4.148X10 47 1246 



0.179 8.418X10 45 394 

0.230 2.893X10 46 569 

0.281 6.581 xlO 46 725 

0.331 1.204X10 47 866 

0.382 1.930X10 47 996 

0.433 2.827X10 47 1114 



1500 




4.0X10 10 6.0X10 10 8.0X10 10 1.0X10 11 1.2X10 11 

r (cm) 

Fig. 9.— RSG25 shocks at the limit of the CASTRO simulated 
domain for six different choices of TOV limit, full neutrino model. 
The choice of TOV limit has a stronger effect on the final shock 
strength in this model than it does in RSG15. 

before reaching the hydrogen envelope. 

Additional effects come into play when using the full 
neutrino model based on the speed at which the core 
accretes. RSG25 has a denser inner structure (higher 



7 



1500 




r (cm) 

Fig. 10.— Shock modeled by CASTRO in RSG15 (TOV = 2.5, 
full neutrino losses) mapped into KEPLER. 

then continued to evolve them. Figures \W\ and QT] show 
the KEPLER results for RSG15, TOV = 2.5; Figure [TO] 
shows the imported shock, and Figure [TT1 shows the final 
velocity of the hydrogen envelope at t = 5 x 10 7 s. 

The shock has decreased significantly in strength by 
the time it reaches the base of the hydrogen envelope; 
however, this envelope is very tenuously bound in both 
RSG15 and RSG25. For each model we tested six choices 
of TOV limit (2.0 - 2.5 M Q , in 0.1 M increments) and 
evaluated the strength of the shock that reached the hy- 
drogen envelope. Using the full neutrino loss model, 
we found in every case tested for RSG15 and in 3 of 
6 tested for RSG25 that the shock produced was larger 
than lxlO 47 ergs, the approximate binding energy of the 
envelope (see Table 3). We can therefore realistically ex- 
pect the envelope to be ejected in these cases. However, 
the highest kinetic energy achieved was only of the or- 
der of 6xl0 47 ergs, and most models fell well below that 
number. The envelope is therefore ejected with a very 
low velocity (50 - 100 km/s). It emits most of its energy 
via hydrogen recombination. Optically this transient has 
a low luminosity ~ 10 39 — 10 40 ergs/s, but maintains this 
luminosity for of order a year. The color temperature of 
the transient is very red, of order 3000 K. An example 
light curve can be seen in Figure [T2l for RSG15, TOV = 
2.5. 

The transients calculated here are obviously much 
fainter and less energetic than standard core-collapse su- 
pernova, but they do bear some similarity to a class of 
recently-observed tr ansients: the "luminous red novae," 
such as V838 Mon (|Munarill2002lh Luminous red novae 
are too bright to be ordinary classical novae, but too faint 
and red to be supernovae. Although V 838 Mon is now 
suspected to be a stellar merger event (|Tvlendall2011| ). 
these two mechanisms have similar end results: a massive 
hydrogen envelope ejected at low energies. Spectroscopic 
observations show however that these phenomena have 
dispersion velocities significantly higher than calculated 
here. The observation of further transients may decide 
this question, or a search for remnants. The shedding 
of a common envelope by a binary merger would leave 
behind a degenerate remnant, but a failed core-collapse 



150 



E 




-o U , , , i , , , i , , , i , , , i 

4 6 8 10 12 

Enclosed Mass (M sun ) 



Fig. 11. — Velocity of the hydrogen envelope at 5x 10 7 s after core 
collapse in RSG15, TOV = 2.5, full neutrino loss model, evolved 
further in KEPLER. 
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Fig. 12.— KEPLER light curve for a transient from RSG15, 
TOV = 2.5. The transient is low luminosity but lasts for around 
a year. 

explosion wo uld leave a black hole. A survey such as that 
proposed by IKochanek et ail (|2008[ ). monitoring red su- 
pergiants for anomalous transients that might signal the 
birth of a black hole, should catch these events. They 
would be visible as a sudden brightening of the "star" 
for of order a year, followed by a gradual but complete 
disappearance. 

In RSG25 a TOV limit of 2.2 M or lower resulted 
in such weak outgoing shocks that they could not be 
accurately followed using KEPLER, and would probably 
be unable to eject the envelope. In situations where the 
envelope is not ejected, there is still the possibility of a 
transient at late times if the envelope is rotating. As it 
falls back into the black hole, the massive envelope may 
create a disk and potentially a long-duration gamma - 
ray transient as described by iWooslev fc Hegen (|2011l) . 
Since the most massive stars are the ones more likely 
to produce black holes quickly, it remains possible to 
produce these long gamma-ray transients. This type of 



transient, while invisible in the optical, could emit low 
levels of gamma rays for months. 

A higher TOV limit in the neutron star EOS will 
increase the probability of these transients occurring. 
Holding the TOV limit constant, the final strength of 
the shock is highest in stars with both smaller initial 
iron core masses (more time spent as a neutron star) 
and smaller carbon-oxygen and helium core masses. We 
might therefore expect the strongest transients to come 
from the lowest-mass red supergiants that fail to form 
CCSNe. Nucleosynthetic constraints place a lower limit 
on the maximum mass star that must exp l ode as a super- 
nova most of the time. iBrown fc Wooslevl (|2013[ ) sets this 
limit at between 20 and 25 M Q . Stars above 20 M Q be- 
come more difficult t o explode, as measured b y their com- 
pactness parameter (jO'Connor fe Qttl 1201 lh . and hence 
more likely to fail. At the same time, in stars above 25 
M Q , it will become increasingly difficult for the shock to 
reach the surface. We may therefore expect the progen- 
itors of these transients, if they do occur, to land in the 
range 20 - 25 M Q . 

For heavier stars that lose their hydrogen envelope and 
die as WR stars, or for stripped progenitors in binaries, a 
shock can form and may reach the surface, depending on 
the size of the remaining helium core. Without a large 
envelope to eject, the transient will be brief but brighter. 

5. CONCLUSIONS 

We have demonstrated that iron core collapse in a mas- 
sive star is capable of producing a faint observable tran- 
sient even if the collapse itself creates no prompt outgo- 
ing shock. The mass lost to neutrinos results, in some 
cases, in unbinding the hydrogen envelope. The amount 
and history of the neutrino mass loss has a strong effect 
on the magnitude of the shock produced, as does the 
structure of the carbon-oxygen and helium cores of the 
progenitor star. In the two red supergiant models tested, 
the shock reached the base of the hydrogen envelope in 
a majority of the models with enough energy to eject it. 
These unusual transients will appear as low-energy, long- 
duration, red events as the ejected envelope emits its en- 
ergy via hydrogen recombination. The ejected envelope 
has a speed on the order of 50 - 100 km/s and maintains 
a luminosity 10 39 — 10 40 ergs/s for approximately a year. 

For a given parametrization of the neutrino losses, the 
transient produced becomes weaker as the TOV mass 
limit is reduced and as the mass of the presupernova 
helium core increases. It therefore remains possible, de- 
pending upon the TOV limit assumed, to fail to eject 
the envelope in more massive stars. If the star in these 
cases has sufficient angular momentum in its outer lay- 
ers, it may in stead produce long gamm a-ray transients as 
described by iWooslev fc Hegerl f|201 lh ; otherwise it will 
disapp ear as an unnova as described bv lKochanek et all 



6. ACKNOWLEDGEMENTS 

We thank Luke Roberts for insight and guidance in 
modeling the neutrino losses from a protoneutron star. 
We thank Ann Almgren and John Bell for providing 
the CASTRO code and Shawfeng Dong, Chris Malone, 
and Haitao Ma for assistance with using it. We thank 
Chris Kochanek and the referee for providing helpful 
comments. This research has been supported by the 



National Science Foundation (AST 0909129), the NASA 
Theory Program (NNX09AK36G), the DOE High En- 
ergy Physics Program (DE-FC02-09ER41438) and two 
UC Lab Fees Research Awards (09-IR-07-117968 and 12- 
LR-237070). 



9 



REFERENCES 



Akmal, A., Pandharipande, V. R., & Ravenhall, D. G. 1998, 
Phys. Rev. C, 58, 1804 

Almgren, A. S. 2010, ApJ, 715, 1221 

Brown, J., Woosley, S. E., 2013. [arXivTl30 2.6973 
Burrows, A., 1988, ApJ, 334, 891 
Burrows, A., & Lattimer, J. M. 1987, ApJ, 318, L63 
Chevalier, R. A. 1993, ApJ, 411, L33 

Demorest, P. B., Pennucci, T., Ransom, S. M., Roberts, M. S. E., 

& Hessels, J. W. T. 2010, Nature, 467, 1081 
Fryer, C. L. 1999, ApJ, 522, 413 

Kochanek, C. S., Beacom, J. F., Kistler, M. D., Prieto, J. L., 
Stanek, K. Z., Thompson, T. A., &; Yiiksel, H. 2008, ApJ, 684, 
1336 

Lattimer, J. M., & Yahil, A. 1989, ApJ, 340, 426 
Lattimer, J. M., & Prakash, M. 2001, ApJ, 550, 426 
Li, X.-D. 2003, ApJ, 596, L199 

MacFadyen, A. I., & Woosley, S. E. 1999, ApJ, 524, 262 
MacFadyen, A. I., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 
410 

McKinney, J. C. 2006, MNRAS, 368, 1561 

Montero, P. J., Janka, H.-T., & Mueller, E. 2011. larXiy-TT0 8.3090 
Munari, U., 2002, A&A 2002, 389, L51 
Nadezhin, D. K. 1980, Ap&SS, 69, 115 



O'Connor, E., & Ott, C. D. 2011, ApJ, 730, 70 
Oppenheimer, J. R., and VolkofT, G. M. 1939, Phys. Rev., 55, 374 
Quataert, E., & Kasen, D. 2011, submitted to MNRAS, MNRAS, 
419, LI 

Quataert, E., & Shiode, J. 2012, arXiv: 1202.5036 

Prakash, M., Bombaci, I., Prakash, M., Ellis, P., Lattimer, J., 

Knorren, R. 1997, Phys. Rept., 280, 1 
Smartt, S. J. 2009, ARA&A, 47, 63 

Tylenda, R., Hajduk, M., Kami?ski, T., Udalski, A., Soszy?ski, I., 
Szyma?ski, M. K., Kubiak, M., Pietrzy?ski, C, Poleski, R., 
Wyrzykowski, ?., Ulaczyk, K., 2011, A&A, 528, A114 

Ugliano, M., Janka, H.-T., Marek, A., k, Arcones, A. 2012, ApJ, 
757, 69 

Weaver, T. A., Zimmerman, G. B., & Woosley, S. E. 1978, ApJ, 

225, 1021 
Woosley, S. E. 1993, ApJ, 405, 273 

Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Reviews of 

Modern Physics, 74, 1015 
Woosley, S. E., & Heger, A. 2007, Phys. Rep., 442, 269 
Woosley, S. E., & Heger, A. 2011, ApJ, 752, 32 



